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We present angle resolved photoemission (ARPES) data from insulating A2CUO2X2 (A=Sr, Ca; 
X=C1, Br), Sr 2 Cu 3 04Cl 2 , and Bi 2 Sr 2 MCu208 (M=Er, Dy) single crystals which illustrate that the 
low energy electronic structure of the half-filled cuprates are independent of the apical atom. By 
performing a photon energy dependent study on CagCuChCk we are able to distinguish between 
features which are intrinsic and those which are a result of the photoemission matrix elements. 
We find that the dispersion is independent of photon energy, while an asymmetry in the electron 
occupation probability across the antiferromagnetic zone boundary is robust to variations in photon 
energy. Finally, we show that the d-wave-like dispersion which exists in the insulator along the 
antiferromagnetic zone boundary does not precisely fit the simple rf-wave functional form near the 
nodal direction. 



I. INTRODUCTION 

One approach to unraveling the mysteries of the high 
T c superconductors is to understand how the electronic 
structure evolves from an antiferromagnetic insulator to 
a superconductor upon doping. The first step is then to 
understand the insulator to metal transition. Although 
Mott qualitatively described how a material, predicted by 
band theory to be a metal, would in fact be an insulator, 1 
it remains unclear as to how the details of the electronic 
structure evolve from a half-filled metal to a Mott insula- 
tor. Questions such as how does the gap open, and what 
happens to the information contained in the occupation 
probability n(k) when strong correlations drive the sys- 
tem insulating need to be answered. The technique of an- 
gle resolved photoemission (ARPES) is a natural choice 
to do this. 

Of course, with every well posed question, there are 
experimental challenges which make it difficult to an- 
swer these questions. La2_ a; Sr a ;Cu04(LSCO) is one 
system which has been successfully grown from half 
filling to heavily overdoped; however, ARPES results 
on these crystals remain somewhat uncertain as illus- 
trated by the extremely broad spectra of La2Cu04. 2 
INcb-xCe^CuO^NCCO) is another system which spans 
a similar doping range, and has attempted to answer 
the above questions; 3,4 however, it is not clear whether 
the results can be generalized to include the hole doped 
cuprates. YBa2Cu306+5(YBCO), which can cover the 



range of half-filling to optimal doping, might appear 
as the next suitable candidate, but the presence of 
a surface state and the existence of one dimensional 
chains has clouded the interpretation of its bulk elec- 
tronic structure. 5 The Bi 2 Sr 2 CaCu20 8+5 (Bi2212) sys- 
tem, with its extremely good cleavage plane, has been 
ideal for ARPES. For precisely this reason, the majority 
of ARPES data on the high T c 's to date have come from 
Bi2212. Unfortunately, high quality Bi2212 crystals at 
very low dopings have not been achieved. On the other 
hand, A2CUO2CI2 (A=Sr,Ca) also cleaves extremely well 
and gives high quality ARPES data comparable to that 
seen in optimally doped Bi2212, but in this case, single 
crystals have only been available at half filling. As a 
result, there exists no perfect system to study the elec- 
tronic structure from the antiferromagnetic insulator to 
the heavily over doped metal in the hole doped cuprates. 

The solution has been to make the reasonable as- 
sumption that the low energy physics of the Bi2212 and 
A2Cu02Cl2(ACOC) systems are identical due to the fact 
that the Cu02 planes are common to both structures. 
Thus the entire doping range can be studied. How- 
ever, these two systems also have several differences. 
Sr 2 Cu0 2 Cl 2 (SCOC) and Ca2Cu0 2 Cl 2 (CCOC) have a 
Cu-O-Cu distance of 3.97A and 3.87A respectively, 6 com- 
pared to 3. 83 A for Bi2212, 7 and they also do not possess 
the orthorhombic distortion and superstructure effects 
which plague Bi2212. Finally, the most striking differ- 
ence in ACOC is that the apical oxygen has been replaced 
with a chlorine atom. It is important to test whether or 
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not these differences can have an effect on the low energy 
electronic structure. 

In this paper we will validate the assumption that low 
energy ARPES data on the oxyhalides are indeed rep- 
resentative of photoemission from a generic, half-filled 
Cu0 2 plane, and thus may be reasonably compared with 
ARPES data on hole doped Bi2212. This will be done 
by showing that replacing the apical chlorine with the 
larger, and less electronegative bromine has no effect on 
the low energy electronic structure. Furthermore, it will 
also be shown that heavily undcrdoped Bi2212 near half 
filling, despite having relatively poor spectral quality, 
is qualitatively consistent with the results on the oxy- 
halides. Finally, ARPES on Sr 2 Cu 3 04Cl2, which con- 
tains an additional copper atom in every other Cu0 2 
plaquette, demonstrates that the lowest lying excitations 
attributed to a Zhang-Rice singlet 8 is surprisingly un- 
affected by even a seemingly large modification of the 
Cu0 2 plane. 

Having shown that the ARPES results of the half filled 
Mott insulating cuprate are not system dependent, we 
will turn our attention in the latter half to extracting 
the physics of the half filled Cu0 2 plane contained in 
the single particle spectral function, A(k, uS). Because 
matrix elements modulate the measured photoemission 
intensity, it is important to distinguish which features of 
the data are intrinsic and which are extrinsic. However, 
photon energy (E 7 ) dependent studies designed specifi- 
cally to test the influence of the matrix element on both 
the dispersion E(k), 9 ' 10 and the intensity n(k), 11 have 
had conflicting reports on the magnitude of variations 
caused by matrix elements. So to test whether or not the 
dispersion and the k dependence of the spectral weight 
are impressive manifestations of matrix clement modula- 
tions, we have performed ARPES n(k) mappings over the 
entire Brillouin zone for 5 different photon energies, and 
examined the (0,0) to (ir, n) cut for 13 photon energies. 
We find that the dispersion is independent of photon en- 
ergy as one might expect, and that with few exceptions, 
the remnant Fermi surface 12 is robust despite observing 
strong variations in spectral weight caused by matrix el- 
ements. 

With an understanding of the matrix element we can 
finally turn our attention to the physics. Specifically, 
we will focus on the <i-wave-likc dispersion found in the 
insulator. 12 This observation allows for a natural connec- 
tion between the d-wave form of the high energy pseu- 
dogap seen in underdoped Bi2212 and the dispersion of 
the insulator as first suggested by Laughlin. 13 This con- 
nection is particularly intriguing in light of the fact that 
the high energy pseudogap and low energy pseudogap 
appear to be correlated. 14 The latter of which is directly 
related to the superconducting gap, thus linking antifer- 
romagnetism which is responsible for the dispersion of 
the insulator to d-wave superconductivity. 

The original ARPES dispersion data on CCOC left 
some ambiguity as to whether or not the detailed dis- 
persion near the node exactly fit the d-wave functional 



form. A linear dispersion at the node of the form 
£(7r/2,7r/2) - E(k) oc ||k - (7r/2,7r/2)|| is highly non- 
trivial and several theories which attempt to connect the 
insulator to the superconductor predict precisely such a 
discontinuity in the derivative of the dispersion at the 
node 13,15,16 . However, the t-J model with next near- 
est neighbor hopping terms t' and t", which correctly 
describes gross features of the dispersion, has a func- 
tional form of cos2k x a + cos2k y a which is analytic at 
k x =k y . To investigate this issue we performed ultra high 
resolution ARPES experiments along the antiferromag- 
netic Brillouin zone(AFBZ) to determine the exact na- 
ture of the dispersion near k x —k y . From this data we 
find that the dispersion near the node is non-linear and 
thus can not be fit by the simple d-wave functional form 
of \cosk x a — cosk y a\. 

The paper is organized as follows. Section II presents 
the experimental details. Sections III through VI present 
photoemission results from the various insulating parent 
cuprates illustrating that the low energy spectra is inde- 
pendent of the apical site. Section VII is a photon energy 
dependent investigation of CCOC. Section VIII compares 
the dispersion of CCOC with the pure d-wave functional 
form. Section IX discusses the experimental results with 
regards to current theoretical understanding after which 
we conclude. 



II. EXPERIMENTAL 

A 2 Cu0 2 X2 (A=Sr, Ca; X=C1, Br), Sr 2 Cu 3 4 Cl2, and 
Bi 2 Sr 2 MCu208 (M=Er, Dy) single crystals were grown 
by a slow cool flux method. 17,18 ARPES experiments 
were performed at beamlincs V-3 and V-4 of the Stanford 
Synchrotron Radiation Laboratory. The beamline V-4 
system is capable of achieving an energy and angular res- 
olution better than 15meV and 0.25°. However, for the 
majority of this work such high resolution is not neces- 
sary due to the extremely broad features which are being 
studied. Thus, the high resolution is only occasionally 
utilized while the remaining data presented uses a reso- 
lution of <70meV and ±1°. Crystals were oriented prior 
to the experiment by Laue back reflection, and cleaved in 
situ at a base pressure better than 5 x 10~ n torr. The 
photon flux was varied to determine whether or not the 
insulating sample was charging. If it was then the sample 
temperature was raised when possible to eliminate this 
effect. A lineshape which is independent of the photon 
flux is taken as proof that there is no charging despite the 
fact that the samples are insulators. Slight charging is 
observed for the Ca 2 Cu0 2 Br 2 and Sr 2 Cu304Cl 2 samples 
at room temperature, but the results were reproducible 
and the spectra simply shifted to higher binding energy 
with increased flux indicating a uniform potential barrier 
formed due to charging. 

We also note that the minimum binding energy of an 
insulator can vary up to leV from cleave to cleave even 
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in the absence of any charging. The reason for this is 
not known with certainty, but is likely due to different 
pinning sites for different cleaves. For the oxy chlorides 
we typically find the centroid of the minimum binding 
energy feature to lie between 0.5 and 0.8eV below the 
chemical potential which is determined by the Ep of a 
reference gold sample in electrical contact. 

III. Sr 2 Cu0 2 Cl 2 AND Ca 2 Cu0 2 Cl 2 

Figure 1 presents a comparison between SCOC and 
CCOC along the high symmetry directions. The two are 
nearly identical. This is to be expected as Sr and Ca are 
isovalent and lie between the CuC>2 planes. From (0,0) 
to (-7T, tt) they show a feature which emerges from the 
background, disperses towards the chemical potential, 
reaches a maximum at (tt/2, tt/2) and then loses weight 
rapidly as it pushes back to higher binding energy. 19 We 
note that the centroid of the lowest energy excitation at 
(tt/2, 7r/2) still lies well below the chemical potential (off 
scale), indicative of the fact that these crystals are Mott 
insulators. Along the (0,0) to (tt, 0) cut a more intense 
and significantly more asymmetric peak is observed. Un- 
der certain experimental conditions it becomes clear that 
this strong asymmetry is due to the presence of a sec- 
ond feature which lies approximately 600meV below the 
main band. 20 The dispersion of this second feature ap- 
pears similar to that of the lower binding energy fea- 
ture, although the intensity is not. The feature along the 
(n, 0) cut does not show much dispersion, and lies ap- 
proximately 350meV below the maximum in dispersion 
at (tt/2, tt/2). In spectra where the peak positions are 
not as clear we will attempt to make a quantitative com- 
parison by taking the minimum of the second derivative 
of the spectra. We do not present any detailed fits of 
the spectra due to a lack of theoretical understanding of 
the broad spectral lincshapes inherent to the half-filled 
cuprates. We note that the broad spectral lineshapes of 
the half-filled cuprates are an intrinsic properties of the 
spectra, and not a result of the sample charging. As 
stated above, this is confirmed in ACOC by varying the 
photon flux and observing no change in the spectra. 

The one significant difference between SCOC and 
CCOC is the observed spectral weight at (tt, 0) in CCOC 
which is suppressed in SCOC. On going from (tt,0) to 
(tt, tt) it can be seen that this weight vanishes quickly. 
Although it is not clear why the matrix elements would 
favor CCOC over SCOC in the (tt, 0) region, it is this 
difference which facilitated the identification of a rem- 
nant Fermi surface in CCOC which was not observed in 
previous studies on SCOC 12 . 

One of the most significant features of the insulator is 
the dispersion between (tt/2, tt/2) and (tt, 0). There are 
two aspects to notice in the right most panels of figure 1. 
The most important is clearly the difference of roughly 
350mcV in dispersion. This is the magnitude of the d- 
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FIG. 1. EDCs of Sr 2 Cu0 2 Cl 2 (a) and Ca 2 Cu0 2 Cl 2 (b) 
along the high symmetry direction indicated in the insets. 
The final panel indicates the magnitude of the d-wave-like 
dispersion seen in the oxyhalides, and the difference in line- 
shape exhibited as a function of k. The measurement con- 
ditions were T=150K, £ 7 =22.4eV and T=100K, £ 7 =25.2eV 
for Sr 2 Cu0 2 Cl 2 and Ca 2 Cu0 2 Cl 2 , respectively. 

wave-like modulation seen in the insulator. The second is 
the strong difference in lincshape. At (tt/2, tt/2) the spec- 
tra show a clear peak which resembles a quasiparticle like 
peak, albeit with a very large width(^300meV), while at 
(tt, 0), the spectra merges more continuously into the high 
energy background, somewhat resembling a step func- 
tion. These will be the benchmarks to which the other 
compounds will be compared in determining whether or 
not the data from SCOC and CCOC are representative 
of a single half-filled Cu02 plane. 



IV. Ca 2 Cu0 2 Br 2 

The biggest difference between ACOC and other high 
T c cuprates is the presence of chlorine as opposed to 
oxygen in the apical site. To see whether or not the 
apical site has an effect on the electronic structure we 
present in figure 2 ARPES data from Ca2Cu02Br2 where 
bromine has replaced chlorine in the apical site. From 
band theory calculations the binding energy of the apical 
orbital is -3.4eV, -2.8eV, and -2.6eV for Ca 2 Cu0 2 Cl2, 
Ca2Cu02Br 2 , and La2Cu04 respectively, 21 which sug- 
gests that CCOB is a suitable candidate to determine 
the effect of the apical chlorine as compared with an api- 
cal oxygen on the low energy electronic structure. On the 
cut from (0,0) to (tt, tt) a dispersive feature is clearly ob- 
served with a minimum in binding energy at (tt/2, tt/2) 
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FIG. 2. a) EDCs of Ca 2 Cu0 2 Br 2 from (0,0) to (tt,tt) 
taken with 16.5eV photons at T=372K. b) Comparison of 
(ir/2, 7r/2)(solid line) and (tt, 0)(dotted line). The dashed lines 
mark the positions of maximum curvature. 




(7t.IT)(jt,0) (71/2,71/2} 



FIG. 3. EDCs of Bi 2 Sr 2 ErCu 2 Os along the high symmetry 
directions as indicated in the cartoons. Below is a plot of 
the second derivative of the above EDCs, from which one can 
trace out the dispersion. S 7 =22.4eV and T=100K. 



and an overall bandwidth of approximately 300 meV. In 
panel (b) the spectra at (tt/2, tt/2) and (ir, 0) are com- 
pared. These data look almost identical to those from 
CCOC which are shown in figure 1. Specifically, the 
lineshape at (tt/2, tt/2) is quite sharp and well defined 
compared with the broad feature seen at (tt, 0). The en- 
ergy difference of 270 mcV (determined by the minimum 
of the second derivative of the corresponding spectra) is 
also consistent with that seen in CCOC. We also took 
a limited amount of data on Sr 2 Cu02Br 2 (not shown). 
The general features of Sr2Cu02Br2 were consistent with 
that of the other oxyhalides, although significant charg- 
ing was observed. The data on the oxybromides suggests 
that the low energy excitations are relatively independent 
of the apical atoms and the observed spectral function in 
ACOC thus originates from the half filled CuC>2 plane. 
Consequently, the comparison of ARPES data on ACOC 
at half filling with Bi2212 data at finite doping seems 
valid. 



V. Bi 2 Sr 2 ErCu20 8 AND Bi 2 Sr 2 DyCu 2 8 

Of course, as mentioned previously, the ideal scenario 
to examine the doping evolution is to study the same sys- 
tem through the entire doping range. Therefore, we also 
present Er and Dy doped Bi2212 crystals grown near 
half filling. Figures 3 and 4 show the EDCs for these 
samples along the high symmetry directions. The low 
energy features are not nearly as well defined as in the 
oxyhalides; however, a clear shoulder does emerge from 
the background. The poor definition of these features is 
consistent with previous studies on underdoped Bi2212 
which show that the low energy excitation spectra be- 



come smeared out as one proceeds toward half filling in 
this system. 22 , 23 This is precisely the reason why ARPES 
data on the oxychlorides has been so valuable to under- 
standing the problem of a single hole in an antiferromag- 
nct. The reason for the relatively poor spectral quality of 
heavily underdoped Bi2212 is an open question, although 
it is possibly simply an issue of sample quality, such as 
lattice strain. 

Examining the low energy spectra from (0,0) to (tt, tt) 
of Er doped Bi2212 we see that the shoulder is most 
pronounced at (tt/2,tt/2) with a minimum binding en- 
ergy of -0.27 eV. From (0,0) to (tt, 0) the shoulder de- 
velops at higher binding energy with a value of -0.49 eV 
at (tt, 0) and then disappears again as one travels from 
(tt, 0) to (tt,tt). Although there is significant ambiguity 
in identifying the binding energy of these features, the 
fact that the shoulder has a dispersion of roughly 0.22eV 
between (tt, 0) and (tt/2,tt/2) can be clearly seen in the 
final panel. The second derivative of the spectra from 
which the above values are attained are also presented 
in figure 3. One can see that the dispersion found with 
this method reproduces the above description. Results 
on Dy doped Bi2212 shown in figure 4 clearly mimic the 
behavior seen in the Er doped Bi2212 including the 0.22 
eV difference in energy between (tt/2,tt/2) and (tt,0). 
The overall features of the dispersion and the resulting 
energy difference of roughly 220 meV between (tt/2, tt/2) 
and (tt, 0) in Er and Dy doped Bi2212 is qualitatively 
consistent with the dispersion and the 350 meV d-wave- 
like modulation seen in CCOC. We suspect the quanti- 
tative difference between Bi2212 and ACOC may simply 
be due to the poor definition of spectral features of the 
former, but we can not rule out the possibility that it is 
an intrinsic difference between the two samples. Never- 
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FIG. 4. EDCs of Bi 2 Sr2DyCu 2 08 along the high symme- 
try directions as indicated in the cartoons. i5 7 =22.4eV and 
T=75K. The dashes lines in the right panel mark the positions 
of maximum curvature. 



theless, even with the poor spectral definition in heavily 
underdoped Bi2212, the electronic structure is still quali- 
tatively consistent between the oxyhalide and Bi2212 sys- 
tems. This result along with our data on CCOB, justifies 
the long standing assumption that the ARPES data from 
SCOC and CCOC are representative of a half filled CuC-2 
square lattice. 



VI. Sr 2 Cu30 4 Cl2: Cu 3 4 PLANE 

Sr 2 Cu304Cl2 is a particularly remarkable example of 
the apparent robustness of the structure of the lowest 
energy excitations. This system deviates from the other 
cuprates due to an additional Cu atom located in every 
other CuC>2 plaquette resulting in an in-plane stoichiom- 
etry of CU3O4. The resulting crystal structure has a unit 
cell twice as large as in S^CuC^C^ and rotated by 45 
degrees (see figure 5). Thus (tt, 0) of the CuC-2 unit cell is 
now equivalent to (tt,tt) in the CU3O4 system, and sim- 
ilarly (7r/2, 7r/2) for Sr 2 Cu0 2 Cl 2 is equivalent to (tt, 0) 
for Sr 2 Cu304Ci2. However, instead of presenting results 
in the new CU3O4 basis, we will continue to present all 
momentum values in the original CuC-2 basis. This will 
make for a simpler comparison between the two systems. 

A sample valence band spectra of S^CusCUC^ at 
(7r/2, 7r/2) is shown in figure 6. In fact it is very sim- 
ilar to the valence band of Sr2Cu02Cl2 at (7r/2, tt/2). 10 
We will however, continue to focus on the "foot" of 
the valence band which contains the low energy excita- 
tions. The following panel examines the low energy ex- 
citation spectra at two equivalent points: (+tt/2, +tt/2) 
and (+7r/2, —tt/2). Since the polarization of the incident 
light is horizontal, the suppression of spectral weight at 
(+7r/2, — 7r/2) relative to (+7r/2, +tt/2) indicates that the 
wave function at {tt/2, tt/2) is odd with respect to a mir- 
ror plane at 45° to the Cu-0 bond direction. This obser- 
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FIG. 5. Real space(a) and Reciprocal space(b) cartoons 
of the unit cell of CU3O4 (dashed line). The unit cell of the 
Cu02 (solid line) is shown for comparison. 



vation is consistent with the state having a d x 2_ y 2 orbital 
character. 

In figure 7, two perpendicular cuts through (tt/2, tt/2) 
are shown. Note that the lowest energy feature is most 
well defined near (tt/2, tt/2), the dispersion is isotropic 
about its maximum at (tt/2, tt/2), and as can be seen in 
the final panel, (tt, 0) lies approximately 320 meV below 
(tt/2, tt/2), although significant uncertainty exists in the 
peak position. These results, including the polarization 
dependence, are identical to the case of ACOC. 10 They 
are also consistent with earlier reports on the isostruc- 
tural compound Ba 2 Cu304Cl 2 taken at a higher photon 
energy. 24,25 We confirm the observation that the spectral 
features about (7r/2, 7t/2) have a remarkable resemblance 
to those in Sr 2 Cu0 2 Cl 2 . It is astounding that the prop- 
erties of the Zhang-Rice singlet seen in ACOC are almost 
undisturbed by the drastic change to the CuC>2 lattice in 
the case of Sr 2 Cu304Cl2. 
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FIG. 6. a) Valence band spectra of Sr2Cus04Cl2 at 
(+7r/2, +7t/2). b) An expanded view of the low energy spec- 
tra. A comparison of the spectra at (+n/2, +7r/2)(solid line) 
and (+7r/2, — 7r/2)(dotted line) indicates that the wavefunc- 
tion responsible for the low energy excitations is odd with 
respect to a 45° line relative to the Cu-O bond direction. 
£-,=16.55eV and T=293K. 



•5 




10.5 11.0 11.5 12.0 10.5 11.0 11.5 12.0 10.5 11.0 11.5 12.0 

Kinetic Energy (eV) 

FIG. 7. a) and b) EDCs of SraCuaC^Cla along two high 
symmetry directions as indicated in the inset, c) Compari- 
son of the spectra at (it/2, it/2) and (tv,0). _E 7 =16.5eV and 
T=293K. 

VII. E 7 DEPENDENCE ON E(k) AND n(fc) 

Having shown that the low energy spectra of a half- 
filled C11O2 plane are independent of the material system 
being studied, we next attempt to separate the features 
of the data which are representative of the single par- 
ticle spectral function, A(k, u), from those which are a 
result of the matrix elements. Unfortunately, extracting 
the single particle spectral function from ARPES mea- 
surements is difficult due to the fact that the measured 
photocmission intensity under the sudden approximation 
is a product of the occupied single particle spectral func- 
tion and the matrix element. In interacting electron sys- 
tems, it is impossible to calculate the matrix elements 
exactly, thus further complicating the ARPES analysis. 
Although we note that symmetry arguments can be very 
powerful in understanding some properties of the matrix 
element, 10 in general it is a function of the experimental 
geometry, photon energy, and the electronic wave func- 
tion. Since its details are not well understood, the objec- 
tive in a given photoemission study must be to focus on 
only those features of the data which are robust against 
variations in the experimental conditions. 

Here we present a photon energy dependence study to 
extract the intrinsic E(fc) and n(fc) structure from the 
raw ARPES data which is modulated by the matrix el- 
ement. On a single cleave we measured EDCs along the 
nodal direction for 13 different photon energies from 14.5 
to 33eV which are shown in figure 8. The morphology 
of the spectra are similar to one another. There are two 
aspects we wish to address. First, the minimum binding 
energy occurs near (ir/2, ir/2) for all the photon energies 
studied, and second, with a few exceptions, the intensity 
profile begins to lose weight before the minimum in bind- 
ing energy is reached (photon energies of 16.5, 32, and 
33eV are the exceptions where the intensity is symmetric 
about (n/2, 7t/2)). We will address the latter point in 



more detail below. 

To check if the dispersion is indeed independent of pho- 
ton energy we plot the peak position versus k for each 
photon energy in figure 9a. The peak positions were 
found by taking the minimum of the second derivative of 
each spectra. One can see that to within our experimen- 
tal limits, which were determined by the reproducibil- 
ity of the dispersion on subsequent scans under identical 
conditions and is given by the width of the blue bar, the 
dispersion is independent of photon energy. This agrees 
with most of the previous reports on SCOC 19 ' 20,26 , in- 
cluding one very detailed, recent study 10 , but contrast 
with the results from Refs. 11 ' 27 , the former of which re- 
port that the minimum binding energy position shifts 
by approximately 10% of the (0,0) to (ir, ir) distance to 
(0.397T, 0.397r) when using 35eV photons, (note that this 
is outside of our momenta error bars) They attribute 
this change to matrix elements where the dependence on 
binding energy varies as a function of photon energy. 

Aside from indicating the expected two dimensional 
nature of the dispersion, our photon energy dependence 
has also clearly resolved the presence of a second com- 
ponent in the low energy electronic structure of the 
half filled insulator. 20 Figure 10 presents EDCs in the 
nodal direction for six photon energies from 16.5eV to 
17.5eV. Aside from the feature typically associated with 
the Zhang-Rice singlet, a second component is observed 
at approximately 600mcV higher binding energy. In the 
nodal direction this feature is most clearly resolved at 
17eV. From this data it is clear that when attempting 
to model the data on the insulator, one can not sim- 
ply treat the high energy spectral weight as a featureless 
background. Although the dispersion of this feature is 
difficult to track it appears to mimic the dispersion of 
the lowest energy feature. The strong influence of matrix 
elements and a lack of understanding of the origin of the 
spectral lineshapes prevents us from extracting detailed 
information on this high energy feature and its evolution 
with photon energy. One possible origin of this feature is 
from string resonances. 28 String resonances occur when 
a hole, created in an antiferromagnetic background, ex- 
periences a confining potential due to the energy cost 
associated with disrupting the antiferromagnetic order 
as it hops away from its original location. The lowest 
energy state for the hole in this potential corresponds 
to the Zhang- Rice singlet, while the first excited state 
is predicted to lie roughly 1.8(J/i) 2 / 3 « 0.5eV higher in 
energy. 29 Although we are not certain of the origin of 
the higher energy feature, the string resonance concept 
provides an intriguing possibility for further study. 

Having examined the dispersion, we now look at the 
more difficult problem of extracting the underlying n(k) 
structure. It is of particular interest to determine 
whether or not the underlying intensity profile of A(k, lu) 
has an asymmetry with respect to the antiferromagnetic 
zone boundary as this would put constraints on the valid 
coupling regime for different models. 28 ' 30 For the major- 
ity of data shown in figure 8 we noticed that the inten- 
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FIG. 8. EDCs from a single cleave of CCOC along F — > (ir, n) for 13 different photon energies indicated in each panel 
respectively. The bold spectra indicates (7r/2,7r/2). The angular separation between the top and bottom spectra in each panel 
is 11.4° with an energy resolution, AE < 50meV. T=200K. 
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FIG. 9. a)i?(k) obtained along the nodal direction from the 
13 photon energies presented in figure 8. E(k) is determined 
by the minimum of the second derivative of the EDCs. b) 
n(k) obtained from the same data by integrating the EDCs 
over an 800meV window. The maxima of each curve were nor- 
malized to each other for display. The shaded bar represents 
the momenta error bars. 



sity of the lowest energy excitation begins to lose weight 
before (tt/2, tt/2). By examining the image plot of a rep- 
resentative set of spectra (as done in figure 11) this ef- 
fect can be seen to be quite dramatic. This has also 
been observed previously by several authors with isolated 
photon energies. 9 ' 19 ' 20 ' 26,27 In figure 9b the n(k) curves 
obtained by integrating the EDCs from -0.5 to 0.3eV rel- 
ative to the valence band maximum for each photon en- 
ergy are shown. In presenting n(k) we implicitly use the 
sudden approximation to extract the momentum distri- 
bution function n(k), from ARPES data via the rela- 
tion n(k) = J A(k,u))f(ui)du> where f(u>) is the Fermi 
function. 31 We note here that the experimental quantity 
measured is not precisely n(k), as the photoemission in- 
tensity is weighted by the matrix element and the integra- 
tion window is limited due to the presence of additional 
bands. The energy integration window was chosen so as 
to minimize the contribution from the second component 
seen at higher binding energy, but we note that the re- 
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FIG. 10. EDCs from a single cleave of CCOC along the nodal direction for _E 7 =16.5eV to _E 7 = 17.5eV as indicated in each 
panel. The bold spectra indicates (7r/2, 7r/2). The angular separation between the top and bottom spectra in each panel is 13° 
with an energy resolution, AE = 70meV. T=293K. 



suits are qualitatively independent of the specific energy 
window chosen. 

Clearly with the exception of 16.5, 32, and 33eV the 
intensity profile peaks well before (tt/2,w/2), and drops 
rapidly as one crosses the antiferromagnetic zone bound- 
ary at (7r/2, 7r/2). The fact that 16.5, 32 and 33eV data 
can differ so dramatically is evidence of what the matrix 
elements can do. If the underlying spectral function were 
symmetric about (7r/2,7r/2) then one would expect that 
the distribution of curves in figure 9b would be evenly dis- 
tributed about (7r/2,7r/2). This is the case for a similar 
study on SCOC, which used five photon energies from 
20 to 246V. 11 They observe the intensity profile to be 
more heavily weighted towards (0,0) at i? 7 =24eV which 
then gradually shifts until at i? 7 =20eV the profile is more 
heavily weighted towards (tt, tt). From this they conclude 
that the true underlying n(k) is symmetric. However, the 
data in figure 9 from 13 different photon energies suggests 
otherwise. 10 curves have an asymmetry more heavily 
weighted towards (0,0), 3 are peaked at (7r/2,7r/2), and 
none are more heavily weighted towards (tt, tt). Combin- 
ing figure 9b with the majority of other single photon 
energy studies on SCOC 9, 19,20,26,27 suggests to us that a 
true asymmetry does exist in the underlying occupation 
probability, n(k). 
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FIG. 11. a) EDCs and b) corresponding image plot of 
CCOC taken along the nodal direction with 23.5eV photons. 
The bold spectra indicates (n/2,n/2), and T=200K. 
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As all of the above studies were only from a single 
cut through the Brillouin zone we now examine a more 
global perspective of the intensity profile. Figure 12b)-f), 
presents relative n(k) patterns of CCOC over a Brillouin 
zone quadrant for five different photon energies while 
panel a) is of an optimally doped Bi2212 sample. Spectra 
are taken at the crosses, values in between were gener- 
ated from a linear interpolation, and except for panel (c), 
the data has been symmeterized about k x =k y and the ge- 
ometry of the experimental setup was identical for each 
data set. The in-plane component of the electric field was 
polarized along the Cu-0 bond direction (This is a 45° ro- 
tation relative to the data presented in figures 8 and 10). 
The first two panels reproduce the original comparison 
of metallic Bi2212 and insulating CCOC from which the 
initial identification of a "remnant Fermi surface" was 
made. 12 For a state with d x 2_ y 2 orbital symmetry the 
suppression of weight as one approaches the line k x —0 
is expected. In Bi2212 the only drop in intensity which 
is not naturally explained by matrix elements is where 
a Fermi surface crossing has occurred. As seen in the 
figure the intensity drop matches the traditional method 
for determining a Fermi surface crossing by following the 
dispersion by eye and is indicated by the dots. This in- 
dicates that the surface of steepest descent in the experi- 
mentally determined n(k) is a valid method of determin- 
ing the Fermi surface. The surprise is that a similar drop 
in intensity is observed in the insulator. Although the 
feature is less well defined here, the striking resemblance 
it bears to the metal suggests that the origin is similar, 
and hence it was qualitatively described as a remnant 
Fermi surface. 12 [Note that the remnant Fermi surface is 
still fully gapped by the large Mott gap.] 

Here we examine the effect of changing the photon en- 
ergy to determine if this feature is robust. From panels 
(b) through (f) one immediately notices that the intensity 
pattern varies tremendously for the five photon energies: 
25.2, 16.5, 29, 32.3, and 41eV. However, the variations 
appear predominantly parallel to the (tt, 0) to (0, tt) di- 
rection, while perpendicular to this there exists relatively 
little variation as we noted above in figure 9. The exact 
shape of the remnant Fermi surface may change, but at 
all photon energies used there is a loss of spectral weight 
as one crosses the approximate antiferromagnetic zone 
boundary from (0, tt) to (7r,0). It may appear that the 
remnant Fermi surface crosses the (0,0) to (tt, 0) cut or 
the (tt, 0) to (it, tt) cut depending on the photon energy 
chosen, but the broadness and the variability due to ma- 
trix elements prevent one from clearly identifying the in- 
tensity profile as either case. However, the data clearly 
suggest that globally, there is an asymmetry in intensity 
stronger towards (0,0) than towards (tt,tt) as one crosses 
the region spanned by the black lines in the Brillouin 
zone which is coincident with the antiferromagnetic zone 
boundary. 

Although the n(k) image plots can provide a wealth 
of information and are extremely good for summariz- 
ing data, it is important to also look at the raw data to 




FIG. 12. (color) Integrated spectral weight. The crosses 
indicate where spectra were taken. Except for (c) the data 
is symmetrized about the k x =k y line. Red is maximum, (a) 
optimally doped Bi2212 at £ 7 =22.4eV. The white hashed re- 
gion indicates the approximate location of the Fermi surface 
determined from n(k). The dots illustrate the position of the 
Fermi surface as determined by the traditional method for an- 
alyzing ARPES data, (b) CCOC shows a striking similarity 
to the metal allowing the identification of the white hashed re- 
gion as a remnant Fermi surface. Comparison of (b) through 
(f) show CCOC taken at photon energies of 25.2, 16.5, 29, 
32.3, and 41eV. The intensity maxima varies between different 
panels, but the loss of intensity as one approximately crosses 
the antiferromagnetic zone boundary is a consistent feature. 
The cumulative boundary of the remnant Fermi surface is 
drawn with black lines on panels (b) through (f). Resolution 
AE < 70meV. 

fully appreciate the information being given by the image 
plots. This is done in Figure 13. Panels (a) and (b) plot 
respectively, the EDCs from (0,0) to (tt, tt) and (tt/2, tt/2) 
to (7r,0), from the data sets used to create the intensity 
maps in figure 12b)-f). We find the spectra are qualita- 
tively similar. This is true even at 41eV where the peak 
is poorly defined throughout the zone. To examine them 
more closely, Figure 14 plots both the dispersion of the 
peak position and n(k) together for all the cuts along the 
antiferromagnetic zone boundary. Similar to the nodal 
direction shown in figure 9a, only slight differences exist 
in dispersion, -E'(k), among the five different photon en- 
ergies. However, the spectral intensity, n(k), varies seem- 
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FIG. 13. ARPES spectra along 2 cuts, a) (0,0) (bottom) to 
(7r,7r)(top) and b) (tt/2, 7r/2)(top) to (n, 0) (bottom), with 5 
different photon energies (16.5, 25.2, 29, 32.3, and 41eV). The 
numbers refer to the percentage distance along each cut. The 
intensities of the features vary, but the dispersion remains the 
same. 

ingly randomly with photon energy. In the extreme case 
between 25.2eV and 29eV the intensity is increasing as 
one approaches (tt/2,tt/2) for the former, and decreasing 
for the latter. One might expect such behavior when the 
underlying spectral function does not possess an asym- 
metry. Meanwhile the asymmetry in spectral weight per- 
pendicular to the antifcrromagnetic zone boundary (see 
figure 9b) appears robust to the strong variations in the 
matrix element. 

In general, the outgoing photoelectrons resulting from 
photons with different wavelengths will necessarily have 
a different overlap with the initial wavefunctions of the 
system, and hence different cross-sections. This is im- 
possible to avoid. Ideally, to eliminate such matrix ele- 
ment effects one would average over all possible photon 
energies, and experimental geometries. Here we have ob- 
served large variations along the antifcrromagnetic zone 
boundary that are clear manifestations of the matrix ele- 
ment, while the asymmetric shift of weight towards (0,0) 
with respect to the (ir, 0) to (0, it) line suggests that this 
is a property of the underlying spectral function. Our 
initial report implied that the underlying n(k) structure 
matched that of the LDA Fermi surface at half filling 
without correlations. 12 It is now clear, that the matrix 
elements are strong enough to make such a precise iden- 
tification very difficult. However, we verify that an un- 
derlying asymmetry exists in the spectral function about 
the antiferromagnetic zone boundary which is robust de- 
spite variations with photon energy. Whether it truly 
lies along the LDA Fermi surface, the antifcrromagnetic 
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FIG. 14. (a) Peak position and (b) n(k) from the spectra 
in figure 13 along the antiferromagnetic zone boundary. The 
data along (n, 0) to (0, tt) has been symmetrized. Again the 
dispersion is independent of photon energy while the intensity 
varies randomly indicating that this modulation in intensity 
is solely due to matrix elements. 

zone boundary, or some other contour is a more difficult 
question to answer. 

VIII. ROUNDED NODAL DISPERSION 

The greatest utility of the remnant Fermi surface con- 
cept as it pertains to the high temperature supercon- 
ductors is that a d-wave-like dispersion was identified in 
the insulator, which could thus provide a natural con- 
nection to the d- wave- like high energy pseudogap. 12 In 
ARPES the pseudogap about (ir, 0) is characterized in 
two ways. The first is by the use of the leading edge 
midpoint, which has the same energy scale as the super- 
conducting gap, and is typically defined in the same way, 
only above T c . However, another energy scale exists in 
the spectra which is an order of magnitude larger, and is 
sometimes refered to as a hump or a high energy pseudo- 
gap. This is the feature which provides a direct connec- 
tion with the dispersion observed in the insulator. 32 The 
remnant Fermi surface identifies a topological contour in 
the Brillouin zone of the insulator, which in this case is 
the antifcrromagnetic zone boundary. By examining the 
dispersion along this contour one can identify a "gap" in 
this system, the dispersion of which along the remnant 
Fermi surface fits the d-wave functional form remarkably 
well. We now extend our earlier study on CCOC to ex- 
amine the precise nature of the dispersion near the nodal 
line k x =k y . 12 A simple d-wave dispersion proportional to 
cos k x a — cos k y a\ would produce a linear dispersion with 
a discontinuous derivative perpendicular to the nodal di- 
rection at k x =k y . This is most easily seen along the an- 
tiferromagnetic zone boundary where the above function 
reduces to | sin(fc x a — 7r/2) | . Such a non-analytic disper- 
sion is non-trivial, and hence the presence or absence of 
such a dispersion is of great significance to theories which 
attempt to unify the antifcrromagnetic insulator with the 
d-wave superconductor. 

In figure 15 we present EDCs taken at 0.6° intervals 
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FIG. 15. (a) EDCs and (b) corresponding image plot ol 
Ca2Cu02Cl2 data taken along the antiferromagnetic zone 
boundary at 0.6° intervals, (it/2, it/2) is indicated by the 
bold EDC and the dotted line in the image plot. Data col- 
lected at T=200K, 25.5eV photons, Resolution AE = 15meV, 
and 0.8° angular resolution 

along the antiferromagnetic zone boundary and through 
(7r/2,7r/2). One observes a smooth round dispersion 
through (7r/2, 7r/2). This is even more evident in the im- 
age plot of the same data. To compare with the d-wave 
functional form we must quantify the dispersion seen in 
the raw spectra. The low energy features seen in the in- 
sulator are inherently very broad with a full width at half 
maximum at (7r/2, 7t/2) of 300meV, which is roughly the 
total dispersion seen in this material. Due to this, and 
the fact that the higher energy spectral weight is also 
dispersive and of unknown origin, the significance of any 
particular fit to the data is questionable. For this rea- 
son we have chosen to quantify the dispersion using three 
methods: the peak maximum, the location of maximum 
curvature, and the leading edge midpoint. We compare 
these quantities against Icosfc^a — cos k y a\ in figure 16. 
The straight line represents the simple ri-wave scenario. 
Near the node it is clear that the dispersion is rounded 
in each case. 

One plausible explanation for the flatness of the dis- 
persion near (tt/2,tt/2) could be the dirty d-wave sce- 
nario. Previous ARPES work on undcrdoped Bi2212 
samples also found a flattened dispersion near the node of 
the superconducting and normal state low energy pseudo 
gaps. 33-35 Earlier works suggested this was consistent 
with a dirty ri-wave scenario since a finite density of states 
with zero excitation energy would result in a flattening 
of the node region. However, Mesot et al. argued that 
flattening due to impurities can be ruled out since pair 
breaking due to impurities should cause the overall mag- 
nitude of the gap to decrease, while they observed the 
opposite. 35 Furthermore, we note that ACOC is very sta- 
ble and most likely very chemically pure as evidenced by 



FIG. 16. Examining the curvature about the nodal line us- 
ing three methods for characterizing the dispersion of the data 
in figure 15 plotted with the original data from ref. 12 . A 
straight line corresponds to a simple d-wave dispersion. Each 
curve was set to zero relative binding energy at the node. 
Note that the d-wave-like gap is in addition to the Mott gap 
which is not shown. The increased scatter as one approaches 
(7T, 0) is indicative of the increased difficulty in tracing the 
dispersion away from (7r/2,7r/2). 

many unsuccessful attempts by many groups to dope the 
oxychlorides under atmospheric conditions. 36 The broad- 
ness of the peaks can not be used as evidence for impu- 
rities since the same feature is almost equally broad in 
Bi2212 samples which contain nearly resolution limited 
quasiparticle peaks in their spectra in the superconduct- 
ing state. 37 Although one can not rule out the possibility 
that the observed dispersion is due to impurities we will, 
for the remainder of this paper, assume that the observed 
rounding is intrinsic to the undoped insulator, and dis- 
cuss the consequences. We defer a full discussion on the 
lincshape of the insulator to another paper. 38 
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FIG. 17. Fits of the data presented in figure 16 to 
A|Scos(20) + (1 — B)cos(6<j})\ where <j) is indicated in the 
figure. B characterizes the deviation from a simple d-wave. 
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In an attempt to quantify the deviation from a sim- 
ple d-wave picture we fit the data to A\B cos(20) + (1 — 
B) cos(6</>)| in figure 17, where (p is defined in the figure, 
cos(20) is approximately | cos k x a — cosk y a\, and cos(6</>) 
is the next allowable harmonic for rf-wave symmetry. An 
increase in the cos(6(/>) weight (smaller B) corresponds to 
a flatter nodal region. Experimentally, we find £?=0.81± 
0.01. This is smaller than any underdoped sample mea- 
sured by Mesot et ai, and thus is consistent with their 
interpretation that the rounding is a result of increasing 
antifcrromagnetic correlations. 35 

Finally we note that the above data was taken at 
S 7 =25.5eV. Identical results were obtained at _E 7 =29eV, 
on a second cleave of CCOC at i? 7 =25.5eV, and on 
SCOC using 22.4eV photons. Improving the angular res- 
olution by a factor of 4 by narrowing the slits of our 
analyzer produced the same results. 



IX. DISCUSSION 
A. System Independence 

By examining A 2 Cu0 2 X 2 , (A=Sr, Ca; X=C1, Br) and 
Bi 2 Sr 2 MCu20 8 (M=Er, Dy) we found that the disper- 
sion of the lowest lying excitations in parent cupratcs 
are independent of the particular system studied. Even 
the results of the CU3O4 system, Sr 2 Cu3C>4Cl 2 , were 
similar to the half-filled Cu0 2 systems. We further 
note that, Nd 2 _ a; Ce a ;Cu04 and La^^Sr^CuC^ near half 
filling have been studied as well. 2 ' 3,39 For Nd 2 CuC>4 a 
shoulder on the edge of the valence band exists and ex- 
hibits a minimum binding energy at (tt/2, 7r/2). 3 Perhaps 
this is the strongest evidence supporting the indepen- 
dence of the low energy electronic structure of half-filled 
cuprates on the apical site as it is completely unoccupied 
in Nd 2 CuC>4. In La 2 CuC>4, a broad shoulder exists at 
(7r/2, 7r/2) at 0.5eV which disperse by roughly 200meV to 
higher binding energy as it approaches (it, 0). 2 ' 39 It is also 
interesting to note that the exchange coupling J, which 
enters into the parameterization of the low energy elec- 
tronic structure, are nearly equivalent in La 2 Cu04(135 
meV) and Sr 2 Cu0 2 Cl 2 (125 meV) as determined from 
inelastic neutron and two-magnon Raman scattering. 40 
All of these observations further support the universal- 
ity to all cuprates of the dispersion of the lowest energy 
feature seen in Sr 2 Cu0 2 Cl 2 and Ca 2 Cu0 2 Cl 2 . 

The role of crystallographic structure on the low en- 
ergy electronic structure and superconductivity is still 
not well established. It was believed that the api- 
cal oxygen or structural distortion of the Cu0 2 planes 
was necessary for superconductivity, but the observa- 
tion of superconductivity in Na-doped Ca 2 Cu0 2 Cl 2 and 
Sr 2 Cu0 2 F 2+d - has since rejected these beliefs. 41,42 In Na- 
doped Ca 2 Cu0 2 Cl 2 no structural distortion was found 
down to 10K while superconductivity was observed as 
high as 26K. Perhaps the apical CI stabilizes a struc- 



ture with an increased Cu-O-Cu distance relative to sys- 
tems with an apical oxygen. Clearly this is not signif- 
icant enough to eliminate superconductivity in CCOC 
with an a-axis lattice constant of 3.87A but may be 
responsible for the fact that SCOC, with a=3.97A has 
not yet been successfully doped into a superconductor. 6 
Although band structure claims that the apical halidc 
ion does not appreciably enter the low energy density of 
states of superconducting oxyhalides, 43 one should note 
that hole-doped superconductivity is observed in com- 
pounds containing an apical anion, while electron-doped 
superconductivity is found in systems missing an apical 
atom. This suggests that the role of the apical atom is 
not completely understood. 

B. E 7 Dependence 

Returning to the photon energy dependence of the 
cuprates, let us first consider the independence of the 
observed dispersion E(k) on the choice of photon en- 
ergy. In principle, there are several reasons which could 
cause E(k) to depend on photon energy which we list be- 
low. We know the cuprates are generally believed to be 
two dimensional electronic systems. In reality the wave- 
functions have some finite z extent, and even if they are 
highly localized will have some finite overlap with neigh- 
boring planes. In fact such a coupling must be present to 
create the observed three dimensional long range mag- 
netic order seen at half filling. Certainly, if the wave 
function overlap were large enough to create a small dis- 
persion as a function of k z a change in dispersion would 
be observed with changing photon energy as this varies 
k z . Alternatively, if the matrix clement had a significant 
dependence in the binding energy of the feature which 
varied with photon energy, a photon energy dependent 
dispersion would also result. Along similar reasoning, if 
there are multiple excitation branches there is no reason 
to expect that the matrix elements dependence on pho- 
ton energy would be the same for each branch. Most 
likely an observed variation with photon energy would 
be caused by a combination of several of these factors. 
However, in figures 9 and 14 we have shown that the dis- 
persion of the insulator is independent of photon energy 
to within our experimental uncertainty. Three dimen- 
sionality may still play a role in causing a small shift in 
dispersion with photon energy, but not surprisingly it is 
a safe assumption to treat the electronic structure of this 
half filled insulator as essentially two dimensional. 

What then is the significance of the observed underly- 
ing intensity modulation along the nodal direction which 
appears robust to variations in photon energy? In the 
context of specific many-body models such as the Hub- 
bard model, there is numerical evidence that a struc- 
ture in n(k) survives even when the on site Coulomb 
U drives the system insulating, even though the dis- 
continuity in n(k) which existed in the metal has been 
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washed out. 28 ' 30 ' This effect is linked to the fact that 
n(k) reflects the underlying Fermi statistics of the elec- 
tronic system, where for the specific case of a two dimen- 
sional square lattice that resembles the CuC>2 planes of 
the cuprates, there is a drop in n(k) across a line that is 
close to the antifcrromagnctic zone boundary 28,30 This 
differs in the generalized t-J models where the structure 
in n(k) is washed out. 30,45 

In the metallic state of optimally doped Bi2212 the 
steepest descent of n(k) gives a Fermi surface consistent 
with traditional ARPES analysis methods despite the 
complication of matrix clement effects as shown in figure 
12. In our earlier paper, 12 the n(k) pattern of the insula- 
tor CCOC was found to be strikingly similar to the n(k) 
pattern seen in Bi2212. This realization, coupled with 
many-body theoretical calculations on various forms of 
the Hubbard model 28 ' 30 , suggests that the insulator pat- 
tern contains information that is related to a Fermi sur- 
face which has been destroyed by strong electron-electron 
interactions thus giving a qualitative concept of a rem- 
nant Fermi surface as the surface of steepest descent in 
n(k). In this paper we have found that the remnant 
Fermi surface acts to emphasize a robust feature which 
we observe in the insulator, although its precise shape 
is uncertain due to the effects of the matrix elements. 
While this may not be a rigorous definition, as the Fermi 
surface is only defined for a metal, this useful concept 
allows a practical connection from the pseudogap seen in 
underdoped cuprates to the properties of the insulator. 13 

Unfortunately, the fact that the antiferromagnetic 
Brillouin zone boundary is quite similar to the under- 
lying Fermi surface means that the roles of magnetism 
and that of the Fermi statistics of the underlying band 
structure in this spectral weight pattern will be uncer- 
tain. Therefore we suggest here a few ideal tests on the 
origin of the n(k) pattern. One idea is to study a system 
whose non-interacting Fermi surface is radically different 
from the magnetic Brillouin zone boundary. Since such 
a system compatible with ARPES may be difficult to 
find, perhaps a more realistic approach will be to study 
a system similar to the cuprates, but where the magnetic 
correlations are much smaller, such that at the measured 
temperature, effects due to magnetic correlations could 
be ruled out. 

We have shown that the loss in intensity as one crosses 
the antiferromagnetic zone boundary is a robust feature 
of the insulator Ca 2 Cu02Cl2 which can not be explained 
solely by matrix clement effects. However, the photon en- 
ergy dependence underscores the qualitative rather than 
quantitative nature of the remnant Fermi surface con- 
cept. We argue that much physics can be learned in spite 
of the adverse effects which matrix elements can have in 
ARPES, as long as care is taken to properly sort out the 
intrinsic versus the extrinsic physics. In particular, the 
resulting connection between the d-wave like dispersion 
and the pseudogap in the underdoped regime is robust. 



C. ci-wave-like Dispersion 

Let us now address the significance of the rounded dis- 
persion at (7r/2,7r/2). Numerical calculations show that 
the t-J model describes the dispersion from (0,0) to (tt, tt) 
well, but incorrectly predicts the energy at (tt/2, tt/2) and 
(tt,Q) to be equal 28 . The addition of t' and t" corrects 
this problem 46 and for realistic values of these param- 
eters, also gives an isotropic dispersion which is scaled 
by a single parameter J. 47 The effect of adding t' and 
t" to calculations is to destabilize the single hole Neel 
state about (tt, 0). The resulting state appears to be 
a spin liquid which may display spin charge separation, 
but this only occurs about (tt, 0). 47 ' 48 Perhaps this ex- 
plains why RVB theories, which assume spin charge sep- 
aration throughout the entire zone, incorrectly predict a 
cusp at (tt/2, tt/2) in the spinon dispersion while still cor- 
rectly describing the overall dispersion on a qualitative 
level. 13 ' 15 

The mean field treatment in a simple SDW picture 
applied to the Hubbard model with up to third near- 
est neighbor hopping fails to describe the data. 47 The 
dispersion along (tt, 0) to (0, tt) is much too great for re- 
alistic parameters. However, the dispersion does con- 
tain a rounded nodal region. If one can argue that one 
should rcnormalize these parameters then it is possible 
to obtain a very good fit to the experimental data as 
done by Nuhez-Regueiro 49 where they also show that 
the magnitude of the overall dispersion decreases with 
doping as has been seen in the high energy pseudogap. 
Alternatively, a diagrammatic expansion of the Hubbard 
model containing t' and t" as well as the QED 3 theory 
also reproduce the observed dispersion. 50,51 Finally, the 
SO(5) theory predicts a | cos k x a — cos k y a\ dispersion for 
the insulator when the superconductor has simple d-wave 
pairing. 16 However, in the projected SO (5), the observed 
flattening of the node in the insulator is consistent with 
a flattened dispersion observed near the node of the su- 
perconducting gap. 52 

D. E fc isotropy of (tt/2, tt/2) 

By combining our investigation of the nodal direction 
with that of the rounding along the antiferromagnetic 
zone boundary we can get a measure of the isotropy of 
the dispersion at (tt/2, tt/2) shown in figure 18. The data 
demonstrates a clear anisotropy in the two directions. 
However, it should be noted that in tracking the disper- 
sion from the EDC's that the total dispersion in both 
directions does not exceed 350meV. It is only due to the 
improved momentum resolution that we get a clear indi- 
cation of the existence of some anisotropy in the disper- 
sion about (tt/2, tt/2). It should be noted that the full 
dispersion in the nodal direction is difficult to extract 
due to a loss of spectral weight at (0,0) which has his- 
torically been attributed to matrix elements. A similar 
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FIG. 18. Anisotropy of the dispersion about (it/2, it/2). 
The nodal direction data is reproduced from figure 9, includ- 
ing all 13 photon energies at once, while the data along the 
antiferromagnetic zone boundary (AFZB) is reproduced from 
figure 17. The dispersion in each case is found by taking the 
minimum of the second derivative of the respective EDCs. 



lack of spectral weight near (n, 0) exists in SCOC, but for 
CCOC the total dispersion along the antiferromagnetic 
zone boundary is easily extracted (see figure 1). The 
resulting 350meV dispersion agrees with the estimates 
determined from SCOC. We note that most models can 
easily account for a small degree of anisotropy by intro- 
ducing additional terms such as next nearest neighbor 
hopping. 



X. CONCLUSIONS 

To summarize, we have studied the Mott-insulating 
cuprates: Sr 2 Cu02Cl2, Ca 2 Cu02Cl2, Ca 2 Cu02Br 2 , 
Sr 2 Cu 3 4 Cl 2 , Bi 2 Sr 2 ErCu 2 8 , and Bi 2 Sr 2 DyCu 2 8 . 
The lowest energy excitations measured by ARPES are 
the same in these compounds as for the prototypical insu- 
lators: Sr 2 Cu0 2 Cl 2 and Ca 2 Cu0 2 Cl 2 . Other studies on 
La 2 -xSr x Cu04 and Nd 2 Cu04 have similar results. 2,3,39 
This shows that the low energy excitations are indeed in- 
dependent of the apical atom, and that the results from 
Sr 2 Cu0 2 Cl 2 and Ca 2 Cu0 2 Cl 2 are truly representative 
of a half-filled Cu0 2 plane. This justifies the connec- 
tion between the high energy pseudogap observed in un- 
derdoped Bi2212 and the d- wave-like modulation seen in 
ACOC. 12 ' 13 Our photon energy dependent study enables 
us to extract the intrinsic spectral function which shows 
that the dispersion is indeed independent of photon en- 
ergy, and that despite strong variations due to the matrix 
elements an asymmetry exists in the n(k) about the line 
from (ir, 0) to (0, ir). We find that the dispersion along 
this contour does not fit a simple <i-wave dispersion, but is 
flattened near (7r/2, 7r/2) consistent with a similar round- 
ing observed in undcrdoped Bi2212 samples. 33-35 Numer- 
ical calculations on the t-t'-t"-J model which match all 
aspects of the insulator dispersion including the rounded 
nodal region, show that the antiferromagnetic correla- 
tions are destroyed near (it, 0), indicating that a spin liq- 



uid picture may be more appropriate for interpreting the 
ARPES data. 47 ' 48 However, generalized t-J models fail 
to reproduce the observed asymmetric intensity profile 
which is reproduced in the Hubbard model. 28,30 Possibly 
the Hubbard model including next nearest neighbor hop- 
ping terms 50 is necessary to capture all the physics of a 
single hole in the half-filled Cu0 2 plane. 
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